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ABSTRACT 

Compact binary white dwarfs (WDs) undergoing orbital decay due to gravitational 
radiation can experience significant tidal heating prior to merger. In these WDs, the 
dominant tidal effect involves the excitation of outgoing gravity waves in the inner 
stellar envelope and the dissipation of these waves in the outer envelope. As the binary 
orbit decays, the WDs are synchronized from outside in (with the envelope synchronized 
first, followed by the core). We examine the deposition of tidal heat in the envelope 
of a Carbon-Oxygen WD and study how such tidal heating affects the structure and 
evolution of the WD. We show that significant tidal heating can occur in the star's 
degenerate hydrogen layer. This layer heats up faster than it cools, triggering runaway 
nuclear fusion. Such "tidal novae" may occur in all WD binaries containing a CO 
WD, at orbital periods between 5 min and 20 min, and precede the final merger by 
10 5 -10 6 years. 

Subject headings: stars: white dwarf - star: binary - hydrodynamics - classical novae 

1. Introduction 

Compact white dwarf (WD) binaries (with orbital periods in the range from minutes to hours) 
are important for several areas of astrophysics. The orbits of these systems decay via the emission 
of gravitational waves, constituting the largest signals for the next generation space-based gravita- 
tional wave interferometers. Systems of sufficiently short orbital period will merge within a Hubble 
time, the result of which may produce a variety of exotic objects, such as helium-rich sdO stars, 
R CrB stars and AM CVn binaries. Most importantly, when the total binary mass is near the 
Chandrasekhar limit, the merged WDs may collapse into a neutron star or explode as a Type la 
supernova (e.g., Webbink 1984; Iben & Tutukov 1984). Recent studies have provided support for 
such "double degenerate" progenitors of SNe la. (e.g., Gilfanov &; Bogdan 2010; Di Stefano 2010; 
Maoz et al. 2010; Li et al. 2011; Bloom et al. 2012; Schaefer & Pagnotta 2012). 

The outcome of a WD binary merger depends on the masses of the WDs and their pre-merger 
conditions (e.g., Segretain et al. 1997; Yoon et al. 2007; Loren-Aguilar et al. 2009; van Kerkwijk 
et al. 2010; Dan et al. 2012; Raskin et al. 2012). Most previous studies of pre-merger binary 
WDs have focused equilibrium tides and considered tidal dissipation in a perameterized way (e.g., 
Mochkovitch & Livio 1989; Iben et al. 1998; Willems et al. 2010; Piro 2011). None of these studies 
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have sought to predict the magnitude and location of tidal heating due to dynamical tides, which 
dominate the tidal responses of the binary WDs. 

In two recent papers (Fuller & Lai 2011, 2012, hereafter paper I and paper II, respectively), we 
presented the first ab initio calculations of dynamical tides in realistic WD models. In paper I, we 
considered resonant excitations of WD g-modes during binary decay and showed that the modes 
reach non-linear amplitudes near the surface of the star. This implies that, rather than exciting 
discrete g-modes, the binary companion will excite a continuous train of gravity waves, which 
propagate outward and dissipate in the outer envelope of the WD. We studied such continuous 
tidally excited waves in paper II. For a canonical Carbon/Oxygen WD (consisting of CO core 
with a He-H envelope), we showed that the outgoing waves are primarily launched at the CO/He 
transition region, and propagate toward the WD surface, where they are likely dissipated through a 
combination of non-linear processes and radiative damping. We computed the energy and angular 
momentum flux carried by the waves in order to predict the orbital and spin evolution of WDs 
in compact binaries. We found that such dynamical tides cause the binary WDs to be nearly 
synchronized prior to merger. Furthermore, the tidal heating rate can be quite large at short 
orbital periods (exceeding tens of solar luminosities just before merger, depending on the system 
parameters), potentially leading to significant observable signatures. 

In this Letter, we show that tidal heating may trigger a thermonuclear runaway hydrogen fusion 
event in a CO WD. The observational consequence of such an event would likely be an outburst 
that resembles a classical nova. We call this new phenomenon a "Tidal Nova" (TN). Unlike all 
other types of novae or supernovae, a TN does not rely on mass accretion or collapse. We present a 
simple two-zone model for the angular momentum evolution of a differentially rotating WD, which 
we use to calculate the radial tidal heating profile within the WD. We then evolve the WD model 
including tidal heating to calculate changes in its temperature, luminosity, and internal structure. 
For a wide range of physically plausible parameters, we demonstrate that tidal heating induces a 
thermonuclear runaway event. Finally, we discuss the observational signatures of such an event, 
and compare our predictions to observations of short-period WD binaries. 

2. Energy and Angular Momentum of Tidally Excited Gravity Waves 

Using the method described in Paper II, we calculate the amplitude of tidally excited gravity 
waves inside a WD. We consider a circular orbit with angular frequency fi. The WD spins at an 
angular frequency J) s , and the spin is aligned with the orbit. In the corotating frame, the frequency 
of the dominant I = m = 2 tidal potential is u> = 2(f2 — f2 s ). For a WD of mass M and radius R 
(and given internal structure) with a companion of mass M', the energy and angular momentum 
fluxes carried by the gravity waves can be written as 



j z (n,u) = T (n)F(u) 
E(n,u) = ^TqF(w), 



(i) 

(2) 
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where 

T^) = — (-), (3) 



a 



with $7 = \J GMt I a 3 the orbital angular frequency (Mt = M + M' is the total mass and a is the 
orbital semi-major axis). 

The dimensionless function F(oj) (similar to the tidal lag angle in the language of equilibrium 
tide theory) determines the magnitude of wave excitation, and is strongly dependent on the internal 
structure of the WD and the tidal frequency oj. In Paper II we have calculated F(oj) for O.6M CO 
WD models of various surface temperatures and slow rotation. We found that F(oj) is an erratic 
function of oj because of the "quasi-resonance cavity" formed by the CO core inside the He/H shell. 
However, because of the strong dependence of F(oj) on oj [the envelope of F(oj) approximately 
scales as oj 5 ], at sufficiently short orbital periods, tidal spin-up combined with orbital decay via 
gravitational radiation ensure that oj ~ const. The orbital period at which this transition occurs 
is P c ~ 40 minutes, depending on the WD masses and temperatures [see Eq. (79) of Paper II]. At 
periods P < P c , the tidal energy transfer rate is 

where / is the moment of inertia of the WD, and tew = I a/a I is the binary inspiral time due to 
gravitational radiation, 

— M^)(^)W 3 

When the outgoing gravity waves damp in the WD envelope and locally deposit their angular 
momentum, some of the wave energy is converted into rotational kinetic energy, while the rest is 
converted to heat. The heating rate is 

E bBBt = E(l-^j. (6) 

If the WD maintains some differential rotation, Q s in the above equation should be the rotation 
rate of the layer in which the waves damp, and heat will also be generated through viscous angular 
momentum transport. 



3. Two Zone Model for Tidal Heat Depostion 

Our calculations indicate that the gravity waves reach non-linear amplitudes and break in 
the outer layers of the WD. The location of wave breaking depends on various parameters (e.g., 
orbital and tidal frequencies), but is always at r > 0.92i? and the exterior mass AM < 10~ 4 M 
(Paper II). Since a small fraction of the stellar mass absorbs the entire angular momentum flux, 
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the outer layer spins up rapidly. If it spins up faster than angular momentum can be transported 
to the core, the outer layer will rotate synchronously with the orbit. Outgoing waves approaching 
the synchronized envelope will be absorbed near corotation and deposit their angular momentum, 
causing the synchronized envelope to move to larger depths (see Goldreich & Nicholson 1989). 

We consider a simple two-zone model for the spin evolution of the WD. In this model, the 
envelope of the star rotates synchronously with the orbit (fi e nv = while the core rotates sub- 
synchronously (r^corc < The envelope and core are coupled, with angular momentum being 
transferred to the core according to a parameterized coupling time, t conp . The angular momentum 
of the core-envelope system evolves according to 

"77 (-^env^env) — >7 2 (i2. Uicorc) 7 (^env ^core)i (7) 
UZ tcoup 

"77 (-^core^core) — 7 (^env ^core)> (8) 
(It ^coup 

where I env = I — I CO re is the moment of inertia of the envelope. Here, J z is the angular momentum 
flux which can be calculated from equation Q . We have assumed that the gravity waves are excited 
in the core and absorbed in the envelop^} Consequently, the angular momentum source term j z is 
only present in the envelope evolution equation, although it is dependent on the tidal frequency in 
the core, u core = 2(Q — Q core ). Using f2 e nv = fi 5 equations and Q can be integrated to find I env 
and ^core cis a function of time or orbital period. The mass AM env of the envelope corresponds to 

Ienv ^ (2/3)AM cnv i? 2 . 

The thickness (or AM env ) of the envelope is dependent on the parameter t C oup- In stably 
stratified stars like WDs, angular momentum can be transported by magnetic fields. In the presence 
of a poloidal field B connecting the core and envelope, t COU p can be estimated from the Alfven wave 
crossing time, 

t .4 = ^Wyr(— j (9) 

for our CO WD model. For WDs without an intrinsic magnetic field, angular momentum may 
be transported via the Tayler-Spruit dynamo (Spruit 2002). To estimate t C oup, we calculate the 
effective viscosity for angular momentum transport via the Tayler-Spruit dynamo, vtSi as outlined 
in Spruit 2002^] We find t^s = J^i 7 " /^Ts)dr ~ 10 3 yr (P/45min) 3 / 2 . Thus we expect the coupling 
time to lie in the range t COU p i$ 10 3 yr for the short orbital periods of interest. 



Figure [T] plots the value of AM cnv as a function of orbital period for our 0.6Mq WD model 

^coup 



with a 0.3Mq companion, using values of i C oup ranging from lyr to 10 3 yr. We begin our calculation 



1 This assumption is valid as long as long as the core-envelope boundary is above the C/He transition layer (with 
an exterior mass AM w 10 _2 Mq), which is the region where the outgoing gravity waves are excited. 

2 For simplicity, we have calculated the viscosity vts without including the effects of composition gradients in the 
WD [see equation (32) in Spruit 2002]. A more realistic estimate of the rotational profile of the WD should take 
composition gradients into account. 
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Fig. 1. — The mass AM env of the synchronized envelope as a function of orbital period for a 0.6Mq 
CO WD model with a O.3M companion. The solid (black) line has t COU p = 1 yr, the dot-dot-dashed 
(green) line t coup = 10 yr, the dot-dashed (orange) line t C oup = 10 2 yr, and the dashed (red) line 
icoup = 10 3 yr. 
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at P or b > lhr and use / e nv,o = and f2 C ore,o = 0, as appropriate at long orbital periods where 
tidal effects are negligible. We see that for the range of t coup considered, AM onv remains small 
(< 1O~ 2 M0) at all orbital periods of interest. Thus, the synchronized envelope most likely does 
not extend down to the C/He transition layer where the gravity waves are excited, justifying our 
assumption that J z is a function of fi C ore- However, the envelope does extend to very large optical 
depths, suggesting that binary WDs may be observed to be synchronized at large orbital periods 
even if their cores are not synchronized. Note that since I env <C I, the core of the star contains 
most of the angular momentum, and its spin evolves in the same manner as discussed in Paper II. 



4. Tidal Heating and Unstable Nuclear Burning 

In the two-zone model discussed in §3, the total tidal heating rate -Ehcat may be calculated 
from equation ([6]) with ft s = tl com , and the tidal heat is deposited entirely at the base of the 
synchronized envelope where AM = AM env . In a real WD, the heat deposition will occur over 
a range of depths that depends on the details of wave breaking and viscous angular momentum 
transport. For simplicity, here we choose to deposit the tidal heat uniformly per unit mass in the 
synchronized envelope. The heating rate per unit mass, £heat 5 is then 

4eat = for AM > AM env (10) 

Eheat = ^|p L for AM < AM env . 

Although the radial dependence of this heating function is unlikely to be realistic, we find that the 
results below are not strongly dependent on the form of e n eat • 

To understand the effect of tidal heating on the WD properties, we evolve WD models using 



the extra heating term calculated via equation (10). We use the one-dimensional stellar evolution 



code MESA (Paxton et al. 2010) to evolve our WD models, starting from an initial orbital period 
of one hour. We present results for a O.6M CO WD model with a ~ 10 -4 Mq hydrogen shell and 
a O.3M companion. 

Figure [2] displays the surface temperature as a function of orbital period for our tidally heated 
WD. For comparison, we also show the temperature of a non-tidally heated WD and the "tidal 
heating temperature" , defined as 



[ cfT.hcat 



1<at ' (11) 



tieat 



At long orbital periods (P > 45 minutes), the tidal heating has little effect on the surface tempera- 
ture of the WD. At shorter periods (P < 30 minutes), the temperature becomes substantially larger 
due to tidal heating. Several of the curves end abruptly due to the ignition of a thermonuclear 
runaway event, at which point we terminate our evolution calculations. 
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Fig. 2. — The surface temperature of the 0.6Mq CO WD model with a 0.3Mq companion as 
a function of orbital period, for initial temperatures of 5000 K (top) and 10 4 K (bottom). The 
solid black lines are calculated with t C oup = 1 y r while the dashed (red) lines are calculated with 
^coup = 10 3 yr. The dotted lines are calculated for a WD with no tidal heating and the same initial 
temperature. The (blue) dot-dashed lines correspond to equation (11). The (red) stars mark the 
points at which tidal novae occur. The asterisks mark the position of the secondary of the 12.75 
minute binary WD system SDSS J065133+284423 (Brown et al. 2011). 
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For small values of i C oup , the tidal heat is deposited at shallow depths and quickly diffuses to 
the surface such that the luminosity of the WD is L ~ Lq + -E n eat> where Lq is the luminosity of 
a non-tidally heated WD. However, for larger values of t coup , most of the tidal heat is deposited 
deeper in the WD where it cannot quickly diffuse outward. This leads to lower surface temperatures, 
although the internal temperature may increase substantially. 

Figure [3] shows the interior temperature profile of our WD at three different orbital periods, 
using i cc .up = 10 3 yr. At long orbital periods, the temperature profile is similar to that of a non- 
tidally heated WD. As the orbital period decreases, the interior heats up, with the local temperature 
maximum at AM ~ AM env . If the base of the hydrogen layer reaches a temperature of ~ 10 7 K, 
hydrogen burning will be ignited. 

In the depicted model, the layer just above the He/H transition (at AM ~ 1CP 4 M0) is 
composed of largely degenerate hydrogen gas. The ignition of fusion in this layer can thus spark 
a thermonuclear runaway. In general, our calculations show that these tidal novae occur only in 
initially cool WDs (T e g < 1.2 x 10 4 K in the absence of tidal heating). They do not occur in hotter 
WDs because the hydrogen is not degenerate and can burn stably. Also, tidal novae require that 
the waves deposit some of the heat near the base of the hydrogen layer, i.e., 1O -5 M < AM cnv < 
10~ 3 M Q . In our two zone model, such heating occurs for coupling times lOyr < t coup ^5 10 4 yr. 
Overall, we find that the tidal novae occur at orbital periods 5 min < P or t, < 20 min, depending on 
the location of heat deposition, initial temperature of the WD, and companion mass. 



5. Discussion 

We have shown that under rather general conditions (see the last paragraph of Section 4), 
tidal dissipation in compact WD binaries can lead to nova outbursts prior to binary merger or 
mass transfer. While we do not attempt to predict the detailed observational signal of a tidal 
nova (TN), we speculate that it may be very similar to a classic nova. However, in contrast to 
classical novae in CVs, a TN would occur in a compact system with no evidence for mass transfer. 
Our results indicate that a TN would precede the beginning of mass transfer or merger by about 
*gw/4 ~ 10 5 — 10 6 yrs [see Eq. provided the conditions outlined in the previous paragraph 
are satisfied. 

In most classical novae, the initial outburst is followed by a period of stable hydrogen burning 
at near the Eddington luminosity, in which the hydrogen shell of the WD inflates to a radius of order 
Rq. However, the ultracompact nature of the WD system involved in a TN (where a ~ Rq/4) 
may preclude such a phase because the stably burning hydrogen shell would inflate beyond the 
WD's Roche lobe. This shell may then accrete on to the companion star or be ejected from the 
system. Therefore, we expect most of the hydrogen to be burned or ejected during in a TN. In the 
absence of mass transfer to supply a fresh hydrogen, recurrent novae would be unlikely. Thus, the 
occurrence rate of these TN may be comparable to that of WD mergers involving a CO WD. 
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Fig. 3. — Temperature profile of the WD (as a function of exterior mass AM) at orbital periods of 
45 minutes (black), 20 minutes (green), and 12 minutes (red). These temperatures are calculated 
for the O.6M0 WD model with an initial surface temperature of T e fj = 5000 K, a 0.3M Q companion 
and i COU p = 10 3 yr. 
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Our theory can be constrained by comparing the prediction of our tidal heating calculations 
to observed compact WD binaries. The 12.75 minute system SDSS JJ065133+284423 provides the 
best opportunity (Brown et al. 2011). This system is composed of a primary with T e g- = 16400 K 
and mass O.25M , and a secondary with T e g- ~ 9000 K and mass O.55M . Comparison with Figure 
[2] indicates that the luminosity of the secondary is likely dominated by tidal heating. Our result for 
a CO WD with an initial temperature of 5000 K and a value of t C oup = 10 3 yr is most consistent 
with the observed temperature of the secondary. These results indicate that a TN may occur in 
this system in the future. 

In principle, tidal heating may change the structure of the WD enough to alter the dynamics 
of gravity wave propagation. However, we find that this is not the case (i.e., no interior convection 
zone forms), with the exception of a thermonuclear runaway event. Our simple two-zone model for 
the WD obviously needs improvement, and we have neglected the effects of mixing induced by the 
breaking gravity waves and viscous angular momentum transport. If the mixing is strong enough 
to smooth out the WD composition gradients, the dynamics of gravity wave excitation and tidal 
heat deposition may be altered. Furthermore, if the surface hydrogen mixes into the WD interior 
where it burns, the surface hydrogen layer will be gradually depleted and a TN will not occur. 
Observations of the ejecta of classical novae indicate substantial enrichment with core elements, 
although the mixing mechanism is not well understood (Truran 2002). These and other aspects of 
TN in compact WD binaries warrant further study. 

Future observations may be able to test whether TN occur and in turn provide information 
about the tidal processes at work in WD binaries. The observation of a nova-like event in a system 
with no evidence for mass transfer would be strong evidence for the existence of TN and for the 
tidal heating mechanism studied in this paper. Measurements of hydrogen surface abundances 
in compact WD systems could also constrain our theory. The observation of a WD with a thick 
hydrogen envelope in a very tight (P < 5 minutes) detached binary would indicate TN do not 
usually occur. If WDs in tight binaries are observed to have little to no hydrogen on their surface, 
this may indicate that TN have stripped the surface hydrogen, or that the hydrogen has been 
destroyed due to efficient mixing processes. Observations of compact binary WDs detected in 
future surveys may provide opportunities to test these theories. 

We thank Bill Paxton, Lars Bildsten, and Eliot Quataert for useful discussions. JF acknowl- 
edges the hospitality (Fall 2011) of the Kavli Institute for Theoretical Physics at UCSB (funded 
by the NSF through Grant ll-AstrollF-0016) where part of the work was carried out. This 
work has been supported in part by NSF grant AST-1008245, NASA grants NNX12AF85G and 
NNX10AP19G. 
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